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We present an x-ray reflectivity study of wetting at the free surface of the binary liquid metal 
gallium-bismuth (Ga-Bi) in the region where the bulk phase separates into Bi-rich and Ga-rich 
liquid phases. The measurements reveal the evolution of the microscopic structure of wetting films 
of the Bi-rich, low-surface-tension phase along different paths in the bulk phase diagram. A balance 
between the surface potential preferring the Bi-rich phase and the gravitational potential which 
favors the Ga-rich phase at the surface pins the interface of the two demixed liquid metallic phases 
close to the free surface. This enables us to resolve it on an Angstrom level and to apply a mean- 
field, square gradient model extended by thermally activated capillary waves as dominant thermal 
fluctuations. The sole free parameter of the gradient model, i.e. the so-called influence parameter, 
K, is determined from our measurements. Relying on a calculation of the liquid/liquid interfacial 
tension that makes it possible to distinguish between intrinsic and capillary wave contributions 
to the interfacial structure we estimate that fluctuations affect the observed short-range, complete 
wetting phenomena only marginally. A critical wetting transition that should be sensitive to thermal 
fluctuations seems to be absent in this binary metallic alloy. 

PACS numbers: 61.25.Mv, 61.30.Hn, 68.10.-m, 61.10.-i 
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I. INTRODUCTION 

The concept of a wetting transition that was intro- 
duced independently by Cahni and Ebner and Saam^ 
in 1977, has stimulated a substantial amount of theo- 
retical and experimental workSiii^iSii. Due to its criti- 
cal character it is not only important for a huge variety 
of technological processes ranging from alloying to the 
flow of liquids, but has the character of an extraordinar- 
ily versatile and universal physical concept, which can 
be used to probe fundamental predictions of statistical 
physics. For example, as shall be depicted in more de- 
tail later on, it can be employed in the case of wetting 
dominated by short-range interactions (SRW) in order 
to test the predictions for the break-down of mean-field 
behavior and the necessary "transition" to a renormal- 
ization group regime, where this surface phenomenon is 
significantly affected by thermal fluctuations. 

A wetting transition occurs for two fluid phases in or 
near equilibrium in contact with a third inert phase, e.g., 
the container wall or the liquid-vapor interface. On ap- 
proaching the coexistence critical point the fluid phase 
that is energetically favored at the interface forms a wet- 
ting film that intrudes between the inert phase and the 
other fluid phase. In general, this surface phenomenon 
is a delicate function of both the macroscopic thermody- 
namics of the bulk phases and the microscopic interac- 
tions. 



Whereas one of the seminal theoretical worke^ on the 
wetting transition gave a microscopic view on this phe- 
nomenon, experimental results at this same level of detail 
were only recently obtained through application of X-ray 
and neutron reflection and diffraction techniquea^*2ii£. 
Moreover, almost all of these experimental studies dealt 
with systems, like methanol-cyclohexane and other or- 
ganic materials, that are dominated by long-range van- 
der-Waals interactions^. The principal exceptions to 
this are the studies of the binary metallic systems 
galhum-lead (Ga-Pb)^^, gallium- thallium (Ga-Tl)-^^, and 
gallium-bismuth (Ga-Bi) ^■^j for which the dominant in- 
teractions are short-range. 

Here, we will present an x-ray reflectivity study of wet- 
ting phenomena that occur at the free surface in the bi- 
nary metallic liquid Ga-Bi when the bulk demixes in two 
liquid phases, i.e. a Bi-rich and a Ga-rich liquid phase. 
The fact that the wetting geometry pins the liquid/liquid 
(1/1) interface that divides the two liquids to the surface 
allows measurements that show the compositional pro- 
file of the wetting film at the Angstrom level. We will 
show results for the structure of the film as it evolves to- 
wards the gravity-thinned film. Through a combination 
of this structural information and the bulk thermody- 
namics of the system, it was possible to extract detailed 
information on the dominant interaction parameters gov- 
erning the surface phenomenology. We will show that a 
square gradient theory that is combined with the effects 
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FIG. 1: The atomic fraction (c)-teinperature(r) bulk phase 
diagram of Ga-Bi. The symbols indicate coexistence lines of 
Predel's phase diagram—. The lines show the phase bound- 
aries calculated from thermodynamic datai". The dashed 
lines represent the metastable extension of the (1/1) coexis- 
tence line below Tm- The points are: C-bulk critical point, M- 
monotectic point, A,B,D,E-points on the experimental path. 
The insets illustrate the surface and bulk phases. In region 
II the wetting film is 50 A thick and the Ga-rich fluid is 5 
mm thick. The bold circless in the insets symbolize the Bi- 
monolayer. 



of thermally excited capillary waves provides a reason- 
able description of that interface. In fact, we will be able 
to extract the sole free parameter of that model, i.e. the 
influence parameter n. This, in turn, will allow us to dis- 
tinguish between intrinsic (mean-field) and fluctuation 
(capillary waves) contributions to the interfacial struc- 
ture. On the basis of this distinction we will estimate 
the influence of fluctuations on the observed SRW at the 
free surface. 



The paper is structured as follows: In the first section, 
we introduce the bulk phase diagram of Ga-Bi and relate 
its topology to the wetting transitions observable at the 
free surface of this binary alloy. X-ray reflectivity mea- 
surements on the wetting films along different paths in 
the bulk phase diagram will be discussed in the second 
section. The third section will focus on the thermody- 
namics and structure of the liquid-liquid interface. In this 
section we will develop a square gradient theory in order 
to model the concentration profile at the liquid/liquid in- 
terface. Finally, in the last section we will provide a more 
general discussion of the rich wetting phenomenology at 
the free surface of this binary liquid metal. 



II. BULK & SURFACE THERMODYNAMICS 

The bulk phase diagram of Ga-Bi, cp. Fig. ^ was 
measured by Predel with differential thermal analysis^-. 
It is dominated by a miscibility gap with consolutc point 
C (critical temperature Tc = 262. 8°C, critical atomic 
fraction of Ga, Ccrit = 0.7) and a monotectic temperature, 
Tm = 222°C. In the part of region I with T < Tm, sohd 
Bi coexists with a Ga-rich liquid. At Tm, the boundary 
between region I and region II, a first order transition 
takes place in the bulk due to the liquidisation of pure 
Bi. For Tm < T <Tc (region II), the bulk separates into 
two immiscible phases, a high density Bi-rich liquid and 
a low density Ga-rich liquid. The heavier Bi-rich phase is 
macroscopically separated from the lighter Ga-rich phase 
due to gravity. In region III, beyond the miscibility gap, 
a homogeneous liquid is found. 

Following the observation of Perepezkoi^ that on cool- 
ing fine Ga-rich droplets are coated by a Bi-rich solid 
phase, Nattland et. al. studied the liquid- vapor interface 
in region II using ellipsometry. They found that a thin 
Bi-rich film intrudes between the vapor and the Ga-rich 
subphase in defiance of gravitjii^, as implied by the illus- 
tration in Fig. Q] This was a clear example of the critical 
point wetting in binary systems that was described by 
Cahn^. On approaching to point C, the Bi-rich phase 
necessarily becomes energetically favored at the free sur- 
face and as a result it forms the wetting film that intrudes 
between the Ga-rich subphase and the surface. In fact the 
situation is slightly more complicated since x-ray studies 
indicate that throughout region III the free surface is 
coated by a monolayer of pure B>li. More recently we 
showed that on heating the binary liquid a thicker wet- 
ting layer of Bi rich liquid forms between the Bi mono- 
layer and the bulk Ga rich liquid. This appears to be 
an unusual example of complete wetting that is pinned 
to the monotectic temperature Tm- This phenomenon 
was first discussed by Dietrich and Schick in order to ex- 
plain an analogous finding in the binary metallic alloy 
Ga-Pbii*i^. The nature of this apparent coincidence of 
a surface transition with the first order transition in the 
bulk at Tm can most easily be illustrated by a transfor- 
mation of the (c, r)-diagram to the appropriate chemical 
potential-temperature (/i, T)-diagram that is depicted in 
Fig. Ola) for Ga-Bi. The axis are temperature, T, the 
difference (/iBi — A'Ga) of the chemical potentials of the 
single species, and their sum (/iei + /iGa)- In this plot, 
the (l/l)-miscibility gap of Fig. Q] which strictly spo- 
ken is a liquid/liquid/ vapor (1/1/v) coexistence boundary, 
transforms into a (l/l/v)-triplc line extending from Af to 
C (solid circles). At M this triple line intersects with 
another triple line, i.e. the solid/liquid/vapor (s/l/v)- 
coexistence line (solid squares), rendering M a tetra point 
of four phase coexistenceSS. At M a solid Bi, a Ga-rich , 
a Ga-poor phase, and the vapor coexist. 

The thermodynamics by which the surface transition 
at M is pinned by the bulk phase transition is obvious if 
one considers the topology of the (/i,T)-plot in the prox- 
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FIG. 2; (a) The chemical potential (/i)-teniperature (T) 
bulk phase diagram of Ga-Bi. The axes are tempera- 
ture, T, the difference (^bi — A'Ga) of the chemical po- 
tentials of the single species, and their sum (pBi+PGa)- 
The solid symbols indicate the following coexistence lines: 
(solid circles) liquid/liquid/vapor (1/1/v) triple line, (solid 
squares) solid/liquid/vapor (s/l/v) triple line, (small trian- 
gles) metastable extension of the liquid/liquid/vapor (1/1/v) 
triple line below Tm, (open circles) experimental path B-D 
probing complete wetting. The points are: C-bulk critical 
point, M-monotectic point, A,B,D,E-points on the experi- 
mental path. To illustrate the 3-dimensional structure of the 
phase diagram, a projection of the phase boundary lines on 
the (/iBi — /iGa,r) plane is drawn in the plot. The symbols 
surrounded by a circle show where the aforementioned char- 
acteristic points fall on this projection. (b)(A/ini,T) phase 
diagram: (A-M) is the (s/l/v), and (M-C) is the (1/1/v) coex- 
istence line. The path B-D is in the single phase region, and 
M and G are the monotectic and critical points. Inset: effec- 
tive wetting layer thickness d on A^M (squares) and B— >D 
(open circles). The solid line is a fit to the A^M d- values. 
(c) The measured d along the experimental path. 



imity of M. The wetting of the free surface by the Bi-rich 
phase as well as the bulk transition are driven by the 
excess free energy, A/i,„, of the Bi-rich phase over that 
of the Ga-rich liquid phase^. This quantity is propor- 
tional to the distance between the (l/l/v)-triple hne and 
any other line leading off (l/l/v)-coexistence, e.g. the 
(s/l/v)-triple hne (A^M) or the hne B^D in Fig. Efa). 
The wetting thermodynamics is displayed in a slightly 
simpler way by the plot A/im vs. T in Fig. Hfb). In 
this figure for T > Tm, the (1/1/v) coexistence line trans- 
forms into a horizontal straight line that extends from M 
to C. For T < Tm the horizontal dashed line indicates the 



metastable (1/1/v) extension of the coexistence and that 
is above the solid-Bi/Ga-rich/vapor (s/l/v) coexistence 
line that goes from M to A. This illustrates the observa- 
tion by Dietrich and Schick^^ that the path A^ M leads 
to coexistence, and thus complete wetting is dictated by 
the topology of the phase diagram. 

A more quantitative understanding of the surface wet- 
ting phenomena can be developed by analyzing the grand 
canonical potential, fis per unit area A of the surface^: 
fis/^ = d A/i-|-7o e"'^''^. Here, d is the wetting film thick- 
ness, ^ is the decay length of a short-range, exponential 
decaying potential, 70 its amplitude and A an arbitrary 
surface area. The quantity A/z comprises all energies 
that are responsible for a shift off true bulk (1/1/v) co- 
existence, i.e. the aforementioned quantity A/i,„. The 
formation of the heavier Bi-rich wetting layer at some 
height, h, above its bulk reservoir costs an extra grav- 
itational energy Afj,g=gApmh where Apm is the mass 
density difference between the two phases. Minimization 
of fis in respect to d then yields the equilibrium wetting 
film thickness of the Bi-rich phase d = ^ln(7o/A/.t). In 
fact the gravitational energy is only significant in compar- 
ison with the other terms for very small values A/i,„ and 
for most of the data shown in Fig. El the gravitational 
term can be neglected, allowing of Afj, = Afim- Thus, 
one expects a logarithmic increase of the wetting film 
thickness upon approaching M from A that is given by 
d = ^ln(7o/A/im), in agreement with the experimental 
finding that was presented in our recent Physical Review 
Letter on this phenomenon - cp. inset in Fig. [2Ib)i^. 
The slight deviations for small values of A/Zm, as well 
as the finite value of d ~ 50A along the coexistence line 
M^B are due to the gravitational term. Since this ap- 
proach of the (l/l/v)-coexistence line ends in M with its 
four phase coexistence, the phenomenon is properly de- 
scribed as tetra point wetting. As demonstrated here, the 
occurrence of this complete wetting phenomenon at the 
surface is an intrinsic feature of the bulk phase diagram. 

In this paper we shall present additional x-ray reflectiv- 
ity measurements that show the evolution of the wetting 
film on approaching coexistence from point D in regime 
III, and hence along a path probing complete wetting 
that is not dictated by the topology intrinsic to the bulk 
phase diagram, but rather by the experimentator's choice 
of the overall atomic fraction of Ga in the sample Cnom, 
i.e. the path B— +D in Fig. ^and Fig. |lla,b). Further- 
more, we will probe the wetting film structure along an 
on- (1/1/v) coexistence path B— * E. 



III. X-RAY REFLECTIVITY MEASUREMENTS 
A. Experiment 

1. Sample Preparation & Sample Environment 

The Ga-Bi alloy was prepared in an inert-gas box using 
> 99.9999% pure metals. A solid Bi pellet was covered 
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FIG. 3: (a) scattering geometry (b) Sketch of the experi- 
mental setup. The arrow close to the mounting indicates the 
variable position of the thermal shield in respect to the sample 
surface. 



by an amount of liquid Ga required for a nominal con- 
centration Cnom = 88 at% Ga. It was then transferred in 
air into an ultrahigh vacuum chamber. A 24-hour period 
of bake-out yielded a pressure of 10"^" torr. The resid- 
ual surface oxide on the liquid's surface was removed by 
sputtering with Ar+ ions. Using thermocouple sensors 
and an active temperature control on both sample pan 
and its adjacent thermal shield a temperature stability 
and uniformity of ±0.05°C was achieved. A sketch of 
the experimental setup, sample environment resp. can 
be found in Fig. ^h). 

A challenge in all x-ray reflectivity measurements of 
liquid metals is the high surface tension of these sys- 
tems, e.g. pure Ga has a surface tension of w700mN/m. 
It often prevents liquid metal from wetting non-reactive 
container walls or substrates. This also leads to large 
curvatures of the liquid surface hampering x-ray reflec- 
tivity measurements^. Remnant oxide layers that exists 
at the liquid/container, liquid/substrate interface resp. 
even enhance this no n- wetting effect; therefore, we re- 
moved these oxide layers from the Mo sample pan by 
sputtering with Ar+ ions, at a sputter current of 25mA 
and a sputter voltage of 2kV, and could reach wetting 
of the Mo crucible by the liquid metal. This resulted 
in small curvatures of the free surface as was judged by 
eye and later on by x-ray reflectivity measurements. The 
resulting flat surfaces facilitated the accumulation of re- 
liable x-ray reflectivity data sets, particularly for small 
incident angles, a. On the other hand, the wetting of the 
container wall by the liquid alloy promoted some spilling 



of the liquid during the sample movements that were nec- 
essary to track the sample position during the reflectiv- 
ity scan. We circumvented this problem by installing a 
ceramic ring that surrounded the sample pan. The ce- 
ramics is not wetted by the liquid metal and, therefore, 
inhibits any liquid flow at the outermost circumference 
of the sample pan. 

Another experimental challenge is related to x-ray re- 
flectivity measurements from liquids in general. The sur- 
face of a liquid is sensitive to any kind of vibrations or 
acoustic noise. In order to obtain a stable reflected signal 
from the liquid surface, we mounted our UHV chamber 
on an active vibration isolationH2.. 



2. X-Ray Reflectivity 

X-ray reflectivity measurements were carried out us- 
ing the liquid surface reflectometer at beamline X22B at 
National Synchrotron Light Source with an x-ray wave- 
length A = 1.54 A. Background and bulk scattering were 
subtracted from the specular signal by displacing the de- 
tector out of the reflection plane by 0.3°. The scattering 
geometry can be found as inset in Fig. |3| The intensity 
i?(gz), reflected from the surface, is measured as a func- 
tion of the normal component q^ — in/X sin(a) of the mo- 
mentum transfer and yields information on the surface- 
normal structure of the electron density p{z) as described 
by the so-called master formuk?-'^. This formula relates 
the Fourier-transformed electron density gradient with 
the experimentally obtained R/Rp. The symbol Rp de- 
notes the Fresnel reflectivity that is expected from an ide- 
ally flat and sharp surface having the electron density of 
the Ga-rich liquid. The standard procedure for determin- 
ing the electron density proflle p{z) from the measured 
reflectivity R{qz) is to construct a simple and physically 
meaningful density model and to flt its Fourier trans- 
form to experimentally obtained data sets of R/R^^. 
We employ a three-box model??, where the upper box 
represents the Gibbs-adsorbed Bi monolayer the second 
represents the Bi-rich wetting film and the lower box rep- 
resents the bulk liquid. The relative electron densities 
of these boxes correspond to pure Bi (top box) and to 
the electron density of the Bi-rich wetting film, p. The 
quantity psub denotes the electron density of the Ga-rich 
subphase. In the simplest approximation the electron 
density profiles of the interfaces between the different 
phases can be described analytically by error- functions. 
The relevant model for this problem includes three error- 
functions (erf)-diffusenesses for the following interfaces: 
vapor/Bi monolayer, Bi monolayer/Bi rich film, Bi rich 
film/Ga rich bulk phase. The first two interfaces de- 
scribing the monolayer feature were remained unchanged 
in the presence of the Bi-rich film, while only the dif- 
fuseness of the Bi rich film/Ga rich bulk phase inter- 
face, (Tobs, was a variable parameter as a function of T. 
The model also included size parameters that describe 
the thickness of the two upper boxes. The box that de- 




FIG. 4: Illustration of the approximation of the analytic 
electron density profile (dashed line) by thin slabs in order to 
apply the Parratt formalism. Here, only the region close to 
the monolayer feature is depicted. 



scribes the bulk sub-phase extends to infinity. During the 
fitting the thickness of the monolayer box was kept con- 
stant while the thickness of the Bi-rich film was allowed 
to vary. 

The use of the master formula tacitly assumes the va- 
lidity of the Born approximation. This assumption holds 
true for wave vectors, q^,, much larger than the value of 
the critical wavevector Qc of the Ga-rich subphase, where 
multi-scattering effects can be neglected^^i. Here, how- 
ever, we will be interested in features corresponding to 
30-60 A thick wetting films, which when translated to 
momentum-space corresponds to features in q^ around 
0.05-0.1 A that are comparable to the value of the crit- 
ical wavevector Qc ~ O.OSA (The qc-values change only 
slightly depending on T, cp. Table P). In order to deal 
with this we resort to the recursive Parratt formalism.^'* 
for (7z-values close to qc- In this formalism one develops 
a 2x2 matrix that relates the amplitudes and phases for 
the incoming and outgoing waves on both sides of a slab 
of arbitrary dielectric constant. For the present problem 
we approximate the above mentioned three-box profile 
by a large number of thinner slabs whose amplitudes are 
chosen to follow the envelope of the analytic profile - cp. 
fig. 01 Typically the number of slabs was of the order of 
300. 



B. Wetting film structures approaching 

(l/l/v)-coexistence 

Here we will present x-ray reflectivity measurements 
from the free surface of Ga-Bi along a path approach- 
ing (1/1/ v) coexistence from the homogeneous part of the 
bulk phase diagram (regime HI), i.e. the path D^B . 
This path, particularly the position of its characteristic 
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FIG. 5: (a) Normalized reflectivities R/Rf for selected points 
D and B corresponding to rD=255.0°C, Tb = 240. 4°C. 
Dashed lines indicate R/Rf = 1. (b) T-dependent nor- 
malized reflectivities R/Rf while approaching coexistence 
on path D — >B along with fits. Dashed lines indicate 
R/Rf ~ l.(c) Electron density profiles p/paub- AH regions 
with p/psub > 1 indicating Bi-enrichment as compared with 
the Ga-rich subphase are gray shaded. 



point B, is solely determined by the overall atomic frac- 
tion of Ga in the sample pan, Cnom- While cooling this 
path intersects the (1/1/v) triple line (cp. Fig. ^ at the 
point B, corresponding to temperature Tb of the misci- 
bility boundary (cp. Fig. P). For T < Tb the heretofore 
homogeneous bulk liquid phase separates into the heav- 
ier Bi rich liquid that settles towards the bottom of the 
pan and the lighter Ga rich liquid on top. The actual 
position in the (c,T)-plane was chosen from considera- 
tion of the vanishing electron density contrast between 
Ga-rich and Bi-rich phase upon approaching C on the 
one hand or the (1/1/v) triple line (cp. Fig. [2Ia)) en- 
forcing the spinodal demixing of the afore homogeneous 
bulk liquid. In our case, its position in the (c,T)-plane 
was chosen after a consideration of the vanishing electron 
density contrast between Ga-rich and Bi-rich phase upon 
approaching C, on the one hand, and a reasonably "long 
path" on (1/1/v) coexistence on the other hand. The 
compromise was a value of about 12% for Cnom, which 
corresponds to Tb = 240.4°C. 

X-ray reflectivity R{qz) was measured at selected tem- 
peratures on path D— >B. In order to avoid artifacts as- 



sociated with the fact that the evolution of the wetting 
film is governed by slow diffusion processes2Si2Si2i equili- 
bration was monitored by continuous taking repeated re- 
flectivity scans following increments of small temperature 
steps of 0.5°C Typically the measured reflectivity fluc- 
tuated wildly for a couple of hours, after which it slowly 
evolved to a stable equilibrium. The fits (lines) to these 
"equilibrium" R/R-p (points) are shown in Fig[S{a,b), 
and the corresponding p{z) profiles - in Fig. EIc). At 
point D [Td = 255°C], typical of region III, R/R^ ex- 
hibits a broad peak at low q^ as well as an increased 
intensity around q^ — O.SA"^. The corresponding elec- 
tron density profile p{z) obtained from the fit indicates 
a thin, inhomogeneous film of increased electron density 
(as compared with psub) close to the surface along with 
a peak right at the surface. It is consistent with a segre- 
gated monolayer of pure Bi along with a thin layer of a 
Bi-enriched phase. As the temperature is decreased to- 
wards B, the peak is gradually shifting to lower g^ and its 
width decreases. This behavior manifests the continuous 
growth in thickness of the wetting layer upon approach- 
ing Tb and is in agreement with the thermodynamic path 
probing complete wetting: A/im — > on path D^B. 

Our experimental data clearly indicate film struc- 
tures dominated by sizeable gradients in the electron 
density, which contrasts with frequently used "homo- 
geneous slab" models, but is in agreement with the- 
oretical calculations. These range from density func- 
tional calculations via square gradient approximations to 
Monte Carlo simulations for wetting transitions at hard 
.^g^^jg2j28j29 Inhomogeneous profiles have also been ob- 
served experimentally in microscopically resolved wetting 
transitions for systems dominated by long-range van-der- 
Waals interactions^*^. Clearly detailed interpretation of 
the inhomogeneity for the GaBi wetting films will re- 
quire either a density functional analysis, or some other 
equivalent approach. Nevertheless, even a simple model 
approximating the wetting layer by a slab of thickness 
d allows a reliable determination of the surface potential 
governing this complete wetting transition. In order to do 
so, effective film thicknesses d have been extracted form 
the p{z) profiles. In Fig. |2Ic) we show plots of these 
d- values versus both T and A/im. The last plot shows 
the expected logarithmic behavior and allows determina- 
tion of the values for the parameters of the short-range 
potential, i.e. 70 = 4QQmN/m and f = h.AA^. More- 
over, agreement between the new (d,^m)-behavior along 
D —>■ B and the data for the path A^M provides an 
experimental prove that both paths have the same ther- 
modynamic character, i.e. they probe complete wetting: 
A^m -^ on path A-^M and D->B. 



C. Wetting film structures: On (l/l/v)-coexistence 

As discussed before, on cooling from the one phase 
region (III) we reach the (l/l/v)-coexistence at B. The 
data for the measured dependence of R/R-p on q^ along 
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FIG. 6: (a) T-dependent normalized reflectivities R/Rp while 
leaving coexistence on path B ^D along with fits. Dashed 
lines indicate R/Rf = 1. With increasing T, each R/Rf 
is shifted by 1.2. (b) Electron density profiles p/psuh- All 
regions with p/psuh > 1 indicating Bi-enrichment as compared 
with the Ga-rich subphase are gray shaded. 



coexistence between points B and M are shown in Fig. 
El^a) . The existence of two peaks at low q^ indicates the 
presence of a fully formed thick film {^ 50 A) . The solid 
lines indicate the best fit results corresponding to the real 
space profiles that are shown in Fig. IHIb). The best fit 
value for maximum density of the thick film p/psub — 
1.20 agrees reasonably well with the value of 1.21 that is 
calculated form the phase diagram at point B. 

These results are also reasonably consistent with what 
is expected for a gravity limited slab of uniform density. 
Using ^ = 5.4A, (To = 400m7V/M and the known material 
constants that make up pm (see tableQ]bclow) the calcu- 
lated value for d = dg = ^ln(o'o/A/ig) = 15.6^ = 85A. In 
view of the fact that this estimate does not take into ac- 
count the excess energy associated with associated with 
concentration gradients across the interfaces some over- 
estimation of dg is not too surprising. Nevertheless, this 
rough calculation does show that the wetting film thick- 
ness is expected to be on a mesoscopic rather than on 
the macroscopic length scale that has been observed for 
similar wetting geometries in systems governed by long- 
range, dispersion forces^. 

Upon cooling from point B to M, the intensity of the 



first peak of R/Rp increases as expected due to the grad- 
ual increase of the electron density contrast p/psuh that 
follows from the increase in the Bi concentration of the 
bulk phase. Similarly, there is a small shift in the po- 
sition of the peak indicating that the thickness varies 
slightly from 53A near B to 5QA near E that follows 
from the T-dependence of the density contrast Apni(T) 
as estimated from the bulk phase A/3ni(T) - cp. Table 
U The on-coexistence path B^E is too far away from C 
to expect more pronounced effects on the wetting layer 
thickness due to a vanishing density contrast or the in- 
crease of the influence of the criticality on the interaction 
potentialsSiSi. 

Moreover the gravitational thinning of the on- 
coexistence wetting film thickness to a length scale com- 
parable to the range of the exponentially decaying short- 
range interactions prevents testing the influence of long- 
range, van-der-Waals like atomic interactions on the wet- 
ting behavior that is present regardless of whether the 
system is insulating or metallio^SiS. 

Finally, we would like to highlight the unique wetting 
geometry encountered here: The subtle balance between 
the surface potential that favors the Bi-rich liquid phase 
at the surface and the gravitational potential that favors 
the Ga-rich liquid phase above the denser Bi-rich phase, 
pins the (1/1) interface between the two coexisting phases 
close to the free surface. It is this property that makes it 
possible to study the structure of these wetting films. For 
example, the width of the interface between that film and 
bulk, cTobs , is the case of the decay of the Kiessig fringes 
at low q^. This will be discussed further in the following 
section. 



IV. LIQUID/LIQUID INTERFACE 

In this section, we will discuss the microscopic struc- 
ture of the (l/l)-interface separating the gravitationally 
thinned film of the Bi-rich phase from the Ga-rich sub- 
phase along the on-coexistence path B— >E. We shall first 
describe a simple square gradient theory for this inter- 
face and then extend it by the consideration of thermally 
excited capillary waves. We will use this theory in or- 
der to extract the (l/l)-interfacial profile and the (1/1)- 
interfacial tension from our measurements of the average 
surface structure. 



A. Square Gradient Theory 

At the (l/l)-interface between the two phases, the con- 
centration varies continuously from the bulk concentra- 
tion of the homogeneous Bi-rich phase, c^, to the bulk 
concentration of the homogeneous Ga-rich phase, c^^. As- 
suming that the variation in concentration within the in- 
terface is gradual in comparison with the intermolecular 
distance, the excess free energy for the inhomogeneous 
region can be expanded in terms of local variables c{r ) 



and \/c{ r )?*. The Gibbs free energy density cost within 
the interface can then be expressed in terms of a combi- 
nation of a local function g{c{r , T) and a power series 
in ve(r)^8i34i35i36: 



G^N 



g{c{-t),T) + ]^n{Vc{-t)f 



dV (1) 



where N is the number of molecules per unit volume. The 
function g{c{ r , T) is the Gibbs free energy density that a 
volume would have in a homogeneous solution. The next 
term is the leading term in the power series expansion. 
The coefficient k referred to as the influence parameter 
characterizes the effect of concentration gradients on the 
free energy. From first principles, it is related to the 
second moment of the Ornstein-Zernike direct correla- 
tion function and can be replaced by the pair potential 
weighted mean square range of intermolecular interac- 
tions of the system^. Equation ^ can be considered 
the Landau- Ginzburg functional for this problem. 

Applying (^) to the one-dimensional composition 
change c(z) across the interface and neglecting terms 
in derivatives higher than the second gives for the total 
Gibbs free Energy G of the system: 



G^NAq 



9{c{z)) 



1 



dc{z) 
dz 



dz (2) 



The variable z represents the direction perpendicular to 
the interface. The symbol Aq denotes an arbitrary sur- 
face area. The interfacial tension 711 is defined as the 
excess energy of this inhomogeneous configuration in re- 
spect to the Gibbs free energy G of the homogeneous 
liquid of one of the coexisting phases. The reference sys- 
tem we have chosen is the homogeneous Bi-rich liquid 
with concentration c^: 
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G{c{z) - G{c^) 



(3) 



Combining the gradient approximation |2Jl with the ther- 
modynamic definition of the interfacial tension (O yields: 



711 = AT 



+ 00 



A»(c(=)) + i.(^ 



dz (4) 



where the grand thermodynamic potential A.g{c) is given 
by 



Agiciz)) = g{c{z)) - gic') 



(5) 



It follows then from that the interfacial tension is 
a functional of the concentration profile c{z) at the in- 
terface. The equilibrium surface tension is obtained by 
minimizing this functional leading to the following Euler- 
Lagrange equation: 



A5(c) 



d^cjz) 
~dl^ 



(6) 



T{°C) 


c' 


c" 


Apm 


p/f 


'sub (ca/c) 


p/psuh(obs) 


qc 


dp 


d 


O"obs 


(^calc 


225.0 
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2.30 




1.24 
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0.41 


0.90 


2.20 




1.24 
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0.0499 


84.9 
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0.43 


0.89 


2.09 




1.22 


1.22 


0.0500 


85.2 


52 


15.0 


14.9 


240.0 


0.45 


0.88 


1.95 




1.21 


1.21 


0.0502 


85.5 


53 


17.5 


17.25 



TABLE I: Material parameters of the coexisting liquid phases as calculated from the bulk phase diagram as well as electron 
density profiles as obtained from our fits to R/Rf for selected temperatures T along the on-coexistence path B^E. The atomic 
fraction of Ga in the coexisting Ga-rich, Bi-rich liquid phase are c\ c' resp. The symbols p / psuh{calc) , p / psuh{obs) denote the 
calculated, observed electron density ratios resp. 



This differential equation along with the boundary con- 
ditions that the concentration has to vary from c^ on the 
one side of the interface to c^^ on the other side deter- 
mines c{z) unambiguously. In fact, the simplicity of © 
allows its direct integration, yielding an expression for 
z{c): 



z{c) = Zq + 



A5(c) 



dc 



(7) 



where zq and cq represent an arbitrary chosen origin and 
composition. It should be noted, that from a renormal- 
ization of the integrand of equation ^ , one can conclude 
that any characteristic length scale, such as the intrinsic 
width of the interfacial profile, tTintrj has to scale as ^/k. 
An attractive feature of this gradient theory is that 
it relates the profile to n and g{c{ r ,T) without regard 
for the theoretical basis by which g{c{ r ,T) is derived. 
In this paper we will use the extended regular solution 
model that will be presented in the appendix to model the 
homogeneous Free Energy of the binary liquid metal, par- 
ticularly the one used to model its miscibility gap. The 
sole quantity required then in order to calculate the (1/1) 
interfacial profile and tension is an expression for the in- 
fluence parameter k, that precedes the square gradient'^''. 
It follows that the influence parameter can be extracted 
from the measurement of the microscopic structure of the 
interface. For example, from an appropriate choice for k, 
solving of ijnji or ([7|l will yield a profile characterized by 
the intrinsic width fTintr that can be compared with CTobs- 



B. Capillary Wave Excitations on the (1/1) 
Interface 



The formalism presented so far relies on a simple mean- 
field picture that neglects any thermal fluctuations on the 
(1/1) interface. By contrast, the real experiment is sensi- 
tive to not only the mean-field diffuseness of the interface, 
but is also sensitive to the fact that thermal fluctuations 
contribute additional broadening to the fluid interface. 
An intuitive, semiphenomenological approach to handle 
the interplay of these two quantities was initiated by Buff, 
Lovett, and Stillingeri^Si^. Essentially, they imagined the 
interface as if it were a membrane in a state of tension 



characterized by the bare interfacial energy 711, as cal- 
culated in the former paragraph; it sustains a spectrum 
of thermally activated capillary waves modes whose av- 
erage energy is determined by equipartition theorem to 
be k^T /2. By integration over the spectrum of capillary 
waves, one finds the average mean square displacement 
of the interface, CcaDi as 



fceT 

47r7n(K) 



In- 



(8) 



Here (/max is the largest capillary wave vector that can 
be sustained by the interface. For bulk liquids this is 
typically of the order of 7r/molecular diameter; however, 
it does not seem realistic to describe excitations with 
wavelengths that are smaller than the intrinsic interfa- 
cial width as surface capillary waves^2ii4. Thus Uintr de- 
scribes the intrinsic interfacial width gmax = Tr/crintr- The 
value of the quantity qmin is slightly more complicated 
since it depends on whether the resolution is high enough 
to detect the long wavelength limit at which an exter- 
nal potential w(z), such as either gravity or the van der 
Waals interaction with the substrate, quenches the capil- 
lary wave spectrum. If the resolution is sufficiently high 

then (7,nin = i^/L^ where L^ = 7^ ^"2 ■ If it is not 
then the value of cr^ap is limited by the resolution limited 
average length scale L^^s over which the observed fluctua- 
tions can be resolved. In this case g,„i,i = A^j-os = Tr/iics 
where Ag^es is the projection of the detector resolution 
on the plane of the interface. In our case gros is signifi- 
cantly larger than ^v therefore, qaiin = 9rcs = 0.04A^^. 

Using the fact that capillary waves obey a gaussian 
statistics which manifests itself as an erf-function type 
profile for the average interfacial roughness^ together 
with the observation that the intrinsic profiles c{z) are 
also well described by erf-function profiles, the two contri- 
butions to the interfacial width can be added in gaussian 
quadrature '^. Thus we write the total observed interfa- 
cial width, (Tcaic, as 



calc 



O'intr(K) 
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ilni^)f 



(9) 



Comparison of © and © reveals that both, cTintr and 
Ccap depend on k. Thus, the determination of k leads 
to the problem of finding the self-consistent value that 
yields a value for CTcaic that agrees with the experimen- 
tally determined roughness (Jobs- 



Intrinsic Profile and (I/l) interfacial tension of 
Ga-Bi 



We performed a self consistent calculation of k for the 
observed (1/1) interface at point E, where <Johs{E) k, 12 A. 
In order to do so, we solved the Euler-Lagrange Equation 
(jHl^ which yields the intrinsic c(z) profile as depicted 
and compared with an erf-profile in Fig. Ela). From that 
we calculated the corresponding erf-diffuseness <Jintr{E) 
and via numerical integration of ||2J) the corresponding 
interfacial tension ^{E) as well as the resultant capillary 
induced roughness and the overall average diffuseness, 
fcap and (Tintr resp. After starting with a reasonable 
value of K, as judged by a scaling analysis of ©, this 
procedure was iterated until crcaic(^) = crobs(^)- The 
value we obtained was Kopt = 5.02 • \Q~^^N'm? /mol. The 
other quantities of interest at E are: crcap(£') = 10.32A , 
a;nt,{E) = 6.35 A and 711 (^) = 3.3miV/m. 

It is interesting to observe that ((HJ implies (J-mtr oc -, 
whereas from (jHl one obtains (Tcap oc ^/K. Hence, the 
function of (Tcaic(«^) exhibits a minimum for Ucaic at 
a characteristic value, Kmin- It is interesting to note, 
that the Kopt, determined from our measurements, cor- 
responds to this peculiar value Kmin meaning that the 
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system selects from all possible CTintr, Ccap combinations, (dashed line): c(z) = 

the one with the minimum overall width for the (1/1)- 

interface. 

Assuming a T-independent infiuence parameter along 
with the available thermochemical data sets that are 
available over a wide range of the phase diagram al- 
lows calculation of the T-dependent interfacial profiles 
and (1/1) interfacial tension from very low up to very 
high T, close to Tc- In Fig. ^h) the intrinsic (1/1) in- 
terfacial profiles are plotted for selected T's versus the 
reduced temperature, t — \T ~ Tc\/Tc (Here T denotes 
the absolute temperature in Kelvin units). The increas- 
ing width of the interface while approaching Tc {t — + 0) 
is well-demonstrated. In Fig. I^the i-dependent interfa- 
cial tension, ^ii{T), is depicted. As expected, it vanishes 
as i -^ 0. 

Comparison of 711 (r) with other measurements would 
be a useful test of our measurements and interpretation. 
Unfortunately, measurements are rather difficult and we 
are not aware of any such measurements for Ga-Bi. For- 
tunately, it is still possible to estimate 711 (T) using other 
available experimental data sets on this alloy. For ex- 
ample, Predeli^ noted that in the vicinity of C the (1/1)- 
coexistence curve of Ga-Bi can be represented by a func- 
tion of the form \(j)c - 0| = K{Tc - Tf with /3 ^ 0.33, 
where is the volume fraction of Ga, Bi resp. and K 
a constant^^. This behavior suggests that the demix- 
ing transition in Ga-Bi belongs to the same universality 
class as binary nonmetallic liquid mixtures, i.e. Ising 
lattice model with dimensionality D=3. Moreover, it al- 
lows for an estimation of the interfacial tension based on 
a scaling relation for these quantities, i.e. the hypothe- 
sis of two-scale-factor universality (TSFU)'*^ should ap- 
ply. Using this concept along with available T-dependent 
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FIG. 7: (a) Calculated intrinsic profile c{z) (-) for T = Te = 
225°C, t — 0.07 in comparison with the erf-function profile 

'""-"'' erf ( ^ " ,^ with 



2+2 -^^ Vv^'^intr(S) 

o"intr(-B) ~ 6.35A. (b) Calculated intrinsic concentration pro- 
files c{z) for selected temperatures T, reduced temperatures 
t resp. as indicated in the figure. 
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FIG. 8: Comparison of temperature dependent, measured in- 
terfacial width a (■) with calculated interfacial widths (lines): 
( — ) , intrinsic width aintr (■ — ), and cap- 
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FIG. 9: Calculated (1/1) interfacial tension 711 of Ga-Bi as 
a function of the reduced temperature t. The dashed line 
represents the TSFU prediction for Ga-Bi*-^. The symbol iM 
indicates the monotectic temperature at fM(Ga-Bi)= 0.08. 



FIG. 10: Calculated (1/1) interfacial tension 711 of Ga-Pb as 
a function of the reduced temperature t. The solid points 
are our calculated values. The dashed line represents a fit 
to the measurements by Merkwitz et. alj--. The symbol tyi 
indicates the monotectic temperature at fM(Ga-Pb)= 0.3. 



measurements of the specific heat, Kreuser and Woer- 
mann extracted an expression for the T-dependent in- 
terfacial tension in Ga-Bi^, 7ts = 7tsO ^^ ■, where /i is 
a critical exponent with universal character for a bulk 
demixing transition, /i = 1.26 and 7tso = 66 ± 15mN/m. 
The corresponding t-dependent 7ts (— ) is plotted in 
Fig. |5| A good agreement with the i-dependence of 
711 calculated by our square gradient theory is found, 
e.g. at iE = 0.07, Te = 225°C: j{E) = 3.31mN/m, 
-fts{E) = 5.7mN/rrM. 



D. (1/1) interfacial tension of Ga-Pb 

Aside from the Ga-Bi system the only other metal- 
lic system for which detailed temperature dependent 
measurements'*® and estimations'*^ of the (1/1) interface 
are reported, is Ga-Pb . Given the close relationship 
between Ga-Bi and Ga-Pb (similar constituents, iden- 
tic phase diagram topology with consolute point C and 
monotectic point M) , we felt encouraged to also apply our 
gradient theory to this binary system. Moreover, we as- 
sumed that the influence parameter, k, as determined by 
the measurements on Ga-Bi, also provides a reasonable 
value for this metallic system. Available thermochemical 
data sets for the miscibility gap of Ga-Pb^'^ allowed us 
to calculate the t-dependent (1/1) interfacial tension. As 
can be seen in fig. ^| where our calculated values for 711 
are shown in comparison with the experimental values, 
an excellent agreement is found. 

Overall, our results suggest that our model of an in- 
trinsic, mean-field interface broadened by capillary waves 
describes the (1/1) interface reasonably well. One must. 



however, bear in mind the approximation implicit in the 
presented analysis. For instance, we have assumed that 
the free surface of the liquid alloy behaves like a rigid 
wall, whereas it will have its own spectrum of capillary 
waves that could be coupled to the waves at the interface. 
However, since the interfacial tension of the free surface is 
about two orders of magnitude larger than the (1/1) inter- 
facial energy in this metallic system, we believe that the 
rigid wall assumption is a reasonable starting point. Fur- 
thermore, the calculation presented above has the short- 
coming of a T-indepcndent influence parameter k. It 
should be mentioned, though, that theory predicts a di- 
vergence of K close to Tc for real fluids, k « (T— Tc) "'°^ 
for T —f Tq. In fact, it can even be shown that this diver- 
gence is necessary to obtain the correct scaling behavior 
of 711 (T) near 



V. WETTING PHENOMENOLOGY AT THE 
FREE SURFACE 

1. Effects of Fluctuations on Short-Range Wetting 

Our analysis in section 1 of the complete wetting tran- 
sition at B was based on a mean-field (MF) model for 
a SRW transition. Such a model accurately predicts the 
critical behavior (i.e. the critical exponents) of systems 
close to a phase transition, provided fluctuations can 
be neglected. For SRW transitions, however, it can be 
shown that the upper critical dimensionality Du = 3^. 
The value D^ is the dimension beyond which MF theory 
can be applied successfully. If the dimension is smaller 
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than Du fluctuations are important and one has to re- 
sort to renormahzation group (RG) methods to describe 
the critical behavior of a phase transition. Because its 
upper critical dimensionality is exactly equal to 3, the 
SRW transition has received a great deal of theoreti- 
cal attention, since it allows one to explore the regime 
where the MF behavior breaks down due to fluctuations, 
and the RG approach becomes applicable. This break- 
down depends on the so-called fluctuation parameter, 
w = /cb • T / {A-K"f\\^^) (where S,b is the bulk correlation 
length, and 7u is the (1/1) interfacial tension) which mea- 
sures the magnitude of the dominant thermal fluctua- 
tions, the thermally induced capillary waves at the fluid, 
unbinding interface, in our case at the (1/1) interface of 
the coexisting Bi- and Ga-rich liquids. 

For a complete wetting transition RG analysis of the 
divergence of the wetting film thickness yields a changed 
prefactor of the logarithmic divergence, dRc ~ Crg (1 + 
w/2) In (l/A/x) ^i^^i^^ . An estimation of w using the bulk 
correlation length, ^^j estimated from the aforementioned 
TSFU prediction, yields for the wetting transitions at M, 
cjm = 0.3 ± 0.2 and for the observed wetting complete 
wetting transition at B, oj-q — 0.4 ± 0.2. Thus, the pref- 
actors of the logarithmic increase of the wetting film with 
A/z is changed by only w 10%, 20% resp. in comparison 
with the mean-field prediction. This is well within our 
experimental error of about 30% for the determination 
of ^. (The large error bar on i^ is due to its extraction 
from a slab thickness analysis whereas the measured wet- 
ting films have significant inhomogeneous structures). A 
clear distinction between RG and MF behavior cannot 
be drawn in our case. 

A critical wetting transition characterized by the wet- 
ting film formation on-(l/l/v)-coexistence, while chang- 
ing T, with a similar value of u) should show more pro- 
nounced deviations from MF behavior. Thus, we shall 
estimate the characteristic temperature Tw of that tran- 
sition in the following. 



2. The Critical Wetting Transition in Ga-Bi 

It follows from our measurements that the critical wet- 
ting transition has to be hidden in the metastable range 
of the (1/1/ v) coexistence line somewhere below Tm- An 
estimation of Tw is possible by a consideration of the 
Free energy at the surface for the wet (Bi-rich wetting 
at surface) and non-wet (Ga-rich phase at the surface) 
situation. In order to do so, we introduce the spreading 
energy, 0(T). It estimates the preference of the Bi-rich 
phase in respect to the Ga-rich phase at the free surface 
at any temperature and is given by the following ener- 
getic contributions. In the wet-situation, the energy is 
given by the sum of the free surface tension of the Bi- 
rich phase and the tension of the (1/1) interface between 
the Bi-rich phase and the Ga-rich subphase. The non- 
wet situation is characterized by the surface energy of 
the bare Ga-rich liquid phase, j{Ga — rich)\^. We define 



Q{T) therefore in the following way: 

e(r) = -f{Ga - rich)i^{T) - 

{j{Bi-rich)i,{T)+-fn{T)) 
= {j{Ga — rich)u{T) — ^{Bi — rich)u{T)) 

= r(r) - 7n(T) (10) 

In equation (|10(l the quantity r(T) denotes the difference 
in the bare liquid- vapor surface tensions of the coexisting 
phases. 

From the definition of Q{T) it follows that establish- 
ment of the wetting film is energetically preferred if > 
(wet-situation), whereas for 8 < it is unfavorable (non- 
wet situation). The condition 8(T) = marks, there- 
fore, the transition temperature Tw, where the system 
switches on-coexistence from a non-wet to a wet-situation 
or vice versa. 

In order to estimate 6(r), in addition to 711 (T), as ex- 
tracted in the previous section we need measurements or 
estimations of 7(Ga — rich)\^{T) and ^{Bi — rich)iy{T). 
Measurements of these quantities are reported in the 
literature^2iSiS, the one by Ayyad and Freyland that 
employed the noninvasive method of capillary wave spec- 
troscopy is the most recent one. From these measure- 
ments, we get a conservative T-independent estimate for 
r of the order of lOOmN/m. Our calculation of 7ii(T), 
particularly its extrapolation towards very low temper- 
ature t — > 1 which yields 'y\\{OK) = 75mN/m, indicates 
that it is always significant smaller than F. Our analysis 
suggests, therefore, that 8(T) > for all T meaning that 
no critical wetting transition occurs on the metastable 
extension of the (1/1/v) coexistence line. 

This remarkable conclusion is, on the one hand, in con- 
trast to the findings for the few examples of binary metal- 
lic wetting systems studied so far: In Ga-Pb and Ga-Tl, 
Wynblatt and Chatain et al. estimated a value for Tw 
significanttly below the corresponding monotectic tem- 
peratures Tm, but still at final values of the order of 
0.3 • Tc{K), 0.2 • Tc{K) for Ga-Pb, Ga-Tl resp. (TdK) 
denotes the critical temperature of bulk demixing on the 
Kelvin temperature scale.) On the other hand, it reflects 
the general finding that in binary metallic systems with 
short-range interactions the critical wetting transition oc- 
curs at significantly lower temperatures than this is the 
case in organic liquid systems, where Tvv is found to lie 
above 0.5 • Tc{K). This tendency seems to be driven to 
its extreme in the case of Ga-Bi: Tw — • rc(^) = ^K . 

A semiquantitative argument for the absence of a crit- 
ical wetting transtion in Ga-Bi as compared to the two 
other metallic systems investigated so far might relate 
to the fact that in comparison with the others the liq- 
uid/liquid coexistence region for GaBi system is signifi- 
cantly shifted towards the Ga-rich, high surface tension 
phase and, in addition the miscibility gab is much nar- 
rower. As a result the change in concentration across the 
1/1 interface is smaller for the GaBi system than for the 
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others. In view of the fact that one can express 



7n = K / dz 



dc{z)' 
dz 



^{Ac) 



(11) 



this suggests that 7;; could be expected to be smaller for 
GaBi than for the others. If this effect is important, and 
if it is not compensated by an accompanying reduction 
in the value of F, then it is possible that Q{T) = T — 
7h > for all temperatures, thereby precluding a wetting 
transition. 



for the excess Bi in the bulk liquid phase below the sur- 
face to diffuse, and precipitate as bulk solid Bi at the 
bottom of the liquid pan, the bulk liquid would simply 
undergo bulk spinodal phase separation. This is con- 
sistent with our observations. Unfortunately it is very 
difficult to absolutely prove that our films are stable and 
this issue can't be resolved from the existing evidence. 
Nevertheless we find it difficult to understand why solid 
Bi would wet the surface for T < Tm and not wet the 
surface as T approaches Tm from above. 



3. Tetra Point Wetting <^ Surface Freezing 



VI. SUMMARY 



Finally, we would like to comment on the fact that the 
observed wetting by the Bi rich liquid at the GaBi mono- 
tectic point is not necessarily the only phenomena that 
can occur. In principle one expects at this tetra point 
a similar or even more diverse wetting phenomenology 
as has been predicted^^ and found^^ in the proximity of 
a triple point for one component systems. For example, 
another possibility is that the free surface could be wet 
by solid Bi. In the bulk both solid Bi and the Bi rich 
liquid are stable for T > Tm', however, for T < Tm the 
free energy of the Bi rich liquid is larger than that of 
solid Bi. Consequently, if the solid Bi surface phase had 
a favorable spreading energy it would certainly wet the 
free surface for T < Tm- In fact, on the basis of optical 
and surface energy measurements Turchanin et al^SiS 
proposed a surface phase diagram for which a surface 
frozen" phase of solid Bi forms at the free surface of the 
GaBi liquid for temperatures below Tm- Unfortunately, 
we have a problem with this. The first may be purely 
semantic, but the idea of a "surface frozen phase" origi- 
nates in the experiments by Earnshaw et al. and Deutsch, 
Gang, Ocko et a l^'^'^^ on the appearance of solid surface 
phases at temperatures above that of bulk freezing for 
alkanes and related compounds. For the GaBi system 
the crystalline phase of bulk Bi is stable at temperatures 
well above those of at which Turchanin et al made their 
observations. We believe that if their observations do 
correspond to thermal equilibrium phenomena, the effect 
would be more appropriately described as wetting of the 
free surface by solid Bi. Unfortunately, the second prob- 
lem we have is that we do not believe that surface wetting 
by solid Bi is favored by the spreading energy. If wetting 
by solid Bi were favored, as implied by Turchanin et al, 
then the wetting layers of Bi rich liquid that are the sub- 
ject of this paper would have to be metastable, existing 
only because of a kinetic barrier to the nucleation of the 
solid wetting film. In fact we have regularly observed that 
on cooling below Tm the original fluid, smooth, highly re- 
flecting surfaces became both rough and rigid. The effect 
is exacerbated by the presence of temperature gradients 
and by rapid cooling. We interpreted this as the forma- 
tion of bulk solid Bi, rather than wetting, since we never 
saw any direct evidence for the formation of films, rather 
than bulk Bi. We argued that if there not enough time 



In this paper we reported x-ray reflectivity measure- 
ments of the temperature dependence of wetting phenom- 
ena that occur at the free surface of the binary metallic 
alloy Ga-Bi when the bulk demixes into two liquid phases, 
i.e. a Bi-rich and a Ga-rich liquid phase. We character- 
ized the temperature dependence of the thickness and 
interfacial profile, with Angstrom resolution, of the Bi- 
rich wetting films that form at the free surface of the 
liquid on approaching the (1/1/v) coexistence triple line, 
i.e. along a path of complete wetting . The results, which 
are characterized by large concentration gradients, agree 
with density functional calculations for such transitions 
at hard walls. On the basis of this it was possible to deter- 
mine the short-range surface potential that favors the Bi- 
rich phase at the free surface and hence is governing the 
wetting phenomenology. The short-range complete wet- 
ting transition turns out to be only marginally affected 
by thermal fluctuations. According to renormalization 
group analysis the critical wetting transition should be 
more sensitive to critical fluctuations; however, accord- 
ing to our analysis the wetting layer of Bi rich liquid 
should be present for all temperatures. This is in con- 
trast to the results for the other binary liquid metal sys- 
tems for which the wetting transitions exists but since 
the metastable 1/1/v line it falls below the liquidus line 
it is hidden. 

The largest thickness of the gravitationally limited Bi 
rich wetting layers that formed at the free surface at 
(l/l/v)-coexistence were typically of the order of « 50A. 
As a result of the fortuitous balance between the surface 
potential that favors the Bi-rich Phase and the gravita- 
tional potential which favors the lighter Ga-rich phase the 
interface dividing the two coexisting phases is sufficiently 
close to the free surface that x-ray reflectivity was sen- 
sitive to its microscopic structure. We interpreted these 
measurements in terms a Landau type of square gradient 
phenomenological theory from which it was possible to 
extract the sole free parameter of that model, i.e. the 
influence parameter k- To the best of our knowledge this 
is the first time that such a parameter has been directly 
determined from measurement. Furthermore, it was pos- 
sible to make a distinction between the intrinsic width 
of the interface and the additional contributions to the 
width from thermal capillary wave features fiuctuations. 
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On making using of available bulk thermodynamic data, 
along with the extracted influence parameter, we were 
able to calculate the experimentally difficult to get (1/1) 
interfacial tension for a wide temperature range. As a 
test of our methods we performed the same calculation of 
the (1/1) interfacial tension for Ga-Pb with its analogous 
miscibility. We assumed that the influence parameter, k, 
as determined from our measurements on the (1/1) inter- 
facial structure in Ga-Bi, was also applicable to Ga-Pb 
and, using the thermochemical data sets available for Ga- 
Pb, we obtained good agreement between our calculated 
values and macroscopic measurements of the (1/1) interfa- 
cial tension. This suggests that the value of the influence 
parameter extracted in this study might provide a rea- 
sonable value for the larger class binary metallic alloys. 
If this proves to be true, then the value determined here 
can be used, along with the surface potential, to predict 
wetting properties of this larger class of systems. 
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APPENDIX: BULK THERMODYNAMICS 

HereS, we focus on the part of the phase diagram dom- 
inated by the miscibility gap and the monotectic point 
M. In a liquid-liquid (1/1) equilibrium of binary systems, 
consisting of two components (Ga and Bi resp.), at at- 
mospheric pressure p, temperature T, the compositions 
ci, cn of two coexisting bulk phases (I and II) are given 
by the thermodynamic conditions: 

ML(r) = f^'LiT) f^UT) = f^g{T) (A.l) 

where fi is the chemical potential. The experimental con- 
ditions of the experiments that we are going to present 
are constant temperature, pressure and a fixed number 
of moles of each component, therefore we have to resort 
to the Gibbs free energy G in order to describe the ther- 
modynamic equilibrium: 

G = U-TS + PV (A.2) 

G ^ n g^ riGa MGa + J^Bi MBi (A. 3) 

Since the total number of particles in mol n = nca+^Bi is 
constant and only two components have to be considered, 
the system can be described with a molar Gibbs Free 
Energy g(c,T) that depends on only on the molar fraction 



of one component: 

G = n(c/iGa + (l-c)/iBi) (A.4) 
^ g(c,T) = c^Ga + (1 -c)^Bi (A. 5) 

With this notation the phase equilibrium conditions 
(equation IjA.ip transform into the following two equa- 
tions: 



|i'- 
|i=- 



oc 
9ic') - 



9W') 



nil 



(A.6) 

(A.7) 



Having a model for the free energy of a binary sytem, 
it is possible, then, to calculate ci and cn of the coex- 
isting phases by solving of equation l|A.7|) . The stan- 
dard approach to describe a miscibility gap in a binary 
demixing system is the regular solution model for the 
Gibbs free energji^. The resulting miscibility gap has 
a symmetric shape and is centered around the consolute 
point Ccrit ~ 0.5 in the (c,T)-plane, quite in contrast to 
the miscibility gap of Ga-Bi with its asymmetric shape, 
centered around Ccrit=0.7 in the (c,T)-plane. Therefore, 
it is necessary to resort to an extended regular solution 
model. Relying on data sets for Ga-Bi from the Calphad 
initiativai^, we use such a model based on Redlich-Kister 
polynoms L^{T)^'^-^^ and express g(c,T) in the following 
way: 



gic,T) =. c-go{Ga)ic,T) + il-c)-go{Bi) 

+R-T[cln{c) + (l-c)/n(l-(4].i 
+ Agniix (c,T) 
5 



The first two terms correspond to the Gibbs energy of a 
mechanical mixture of the constituents; the second term 
corresponds to the entropy of mixing for an ideal solu- 
tion, and the third term, A^xs is the so-called excess 
term. Here, it is represented by an extension of the usual 
regular solution expression, i.e. a sum of Redlich-Kister 
polynomials, L^{T), as listed in Table ITU 

By solving the nonlinear equations (|A.7p applied on 
gic,T) for the temperature range Tm < T < 262°C 
we could calculate the binodal coexistence line as plot- 
ted in Fig. n The agreement of the calculated phase 
boundaries with Predel's measured phase boundaries is 
excellent. Particularly, the measured consolute point C 
{Tc — 262°C, Ccrit = 0.7 is nicely reproduced by the 
calculated critical values Tc — 262. 8C and c = 0.701. 
Furthermore, the knowledge of g{c,T) allows us to ex- 
trapolate the binodal lines below Tm into the region of 
metastable (1/1) coexistence, and, hence, to get informa- 
tion on the energetics of the metastable Ga-rich, Bi-rich 
phase resp. Similar to the procedure presented above, we 
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V 


Redlich-Kister polynomial Li, (T) 





80000-3389+T 


1 


-4868-2.4342-T 


2 


-10375-14.127-r 


3 


-4339.3 


4 


2653-9.41-T 


5 


-2364 



also calculated the phase boundaries below Tm- Here, we 
used data sets for the Gibbs free energy of pure solid Bi, 
which coexists for T < Tm with a Ga-rich liquid. 



TABLE II: Redlich-Kister polynomials used to model the (1/1) 
miscibility gap of Ga-BiiS*Si. 
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